Bose-Einstein condensates at equilibrium inside a pancake-shaped trap: dimensional 

cross-over in the scattering properties 
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Two-dimensionality of the scattering events in a Bose-Einstein condensate introduces a logarith- 
mic dependence on density in the coupling constant entering a mean-field theory of the equilibrium 
density profile, which becomes dominant as the s-wave scattering length gets larger than the conden- 
sate thickness. We analyze quantitatively the role of the form of the coupling constant in determining 
the transverse profile of a condensate confined in a harmonic pancake-shaped trap at zero tempera- 
ture. We trace the regions of experimentally accessible system parameters for which the cross-over 
between different dimensionality behaviors may become observable through in situ imaging of the 
condensed cloud with varying trap anisotropy and scattering length. 

q , PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj 
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■ I. INTRODUCTION 

, 

A system of dilute hard-disk bosons constrained to move inside a two-dimensional (2D) box is a well-known model in 
statistical mechanics [Q. This system shows very different features from its three-dimensional (3D) counterpart. First 
of all its collisional properties are very peculiar, since in two dimensions the T-matrix vanishes at low momentum and 
I | energy |l],|2| , and the resulting boson-boson coupling constant depends on the density of the system through the loga- 
rithm of the diluteness parameter na 2 , where n is the areal density of the disks and a the hard-core diameter. Another 
important difference with respect to 3D systems is that while at zero temperature a fraction of the bosons undergo 
Bose-Einstein condensation, at finite temperature phase fluctuations destroy the long-range order (in agreement with 
the Mermin- Wagner theorem [0]) and the one-body density matrix decays algebraically to zero at large distances. 
The interacting 2D Bose gas is nevertheless a superfluid also at finite temperature and the nature of the transition 
is known to be of the Kosterlitz-Thouless type jj], the critical temperature being Tkt ~ h 2 n/(m In ln(l/na 2 )) in the 
ultra-dilute limit lnln(l/rta 2 ) 3> 1 f§- In a recent experiment a 2D Bose gas has been realized using a film of gaseous 



Hydrogen on a Helium substrate 



A dilute 2D Bose gas inside a harmonic trap has also been considered by several authors. For the ideal gas the 
. presence of the confinement introduces essential changes in the thermodynamic properties relative to the homogeneous 
case, and Bose-Einstein condensation starts at a finite temperature ksT c « huj^N 1 / 2 where u>± is the frequency of the 
2D harmonic confinement and N the number of particles @ . The interacting system is expected to be Bose-condensed 
at low temperature T < T^, but phase fluctuations should become important at a temperature where the phase 
correlation length becomes smaller than the radius of the cloud Q. The precise behavior of the system at these 
temperatures and the nature of the phase transition have not yet been fully clarified (see e.g. Fernandez and Mullin 
C \ H and references therein). 

. In current experiments on Bose-Einstein condensates of alkali atoms in magnetic or optical traps the anisotropy 
of the confinement can be varied to obtain flatter and flatter (quasi-2D) condensates [[To), with the ultimate aim to 
observe the special features of low dimensionality. It is thus important to realize where in parameters space and which 
O ■ physical aspects of 2D systems will become dominant as the anisotropy is increased. The key question in this respect 
is which is the appropriate model for the boson-boson coupling strength in the regime of cross-over between 3D and 
2D and how to insert it in a self-consistent calculation of the equilibrium properties of a Bose-Einstein condensed 
rS [ cloud. 

In this paper we treat this problem at zero temperature. We evaluate the equilibrium density profiles of the 
condensate under different choices for its physical parameters, ranging from a 3D anisotropic system to a strictly 2D 
one. With increasing anisotropy the system first becomes 2D with regard to the confinement - only the lowest axial 
state is occupied - and then also in its collisional properties. We find that these different regimes can be identified by 
observing the size of the cloud in the radial plane, and we characterize the cross-over in terms of the relevant physical 
parameters. 

The paper is organized as follows. In Sec. [o] we introduce the coupling strengths used to describe the condensate 
from the 3D to the 2D regime and the corresponding non-linear Schrodinger equation (NLSE). The results for the 



equilibrium density profile arc given in Sec. |IIIj, and Sec. [V presents our conclusions and final remarks 
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II. GROUND STATE AND SCATTERING PROPERTIES OF A QUASI-2D CONDENSATE 



We consider a dilute Bose-condensed gas at zero temperature under anisotropic (pancake-shaped) harmonic con- 
finement characterized by the frequencies lo j_ and to z = \cu± with A 1. Within a density- functional approach, the 
ground-state properties of the gas are described in terms of the condensate wave function i/>(r) in the {x,y} plane, 
which is to be determined as the minimum of the local-density energy functional 



E[ip] = / d 2 r 



^\^\ 2 + (V ext -nM 2 +e(nM< 
2m 



(1) 



Here n = \ip\ 2 is the particle density, /i is the chemical potential, and e(n) = gn/2 is the ground-state energy per 
particle of a homogeneous Bose gas with short-range interactions in the mean-field approximation. As we shall see, 
the coupling g can depend on the condensate density: g = g(i/j). The Euler equation SE/Sip* — leads to the NLSE 
for the condensate wave function [[lT] 

-^V 2 ^ + V ext ^ + g(iPM 2 *P = (2) 

The familiar Gross-Pitaevskii equation is recovered by taking a constant coupling g. Extensions of this equation 
beyond mean field and beyond the local-density approximation have been proposed for systems at higher density 

We now need to specify the choice of g for the system of present interest. Quite generally the coupling g is obtained 
microscopically from the effective interaction potential T(k,k',P) in the limit of low energy and momenta. The 
effective interaction can in turn be related to the two-body scattering function f(k,k') — J tyAr)v{r)e~ r d^r in 
dimension d = 3 or 2, with ip%(r) being the outgoing wave function for the relative motion of the two particles and 

v(r) being the interparticle potential jl5[ ]. For a 3D system f(k,k') is the scattering amplitude and g is a constant 
determined by the s-wave scattering length a: 

^h 2 

gsD = a. (3) 

m 

On increasing the anisotropy of a 3D condensate, its physical properties will change first due to the modified shape 



of the confinement |16|, but then also due to the modified scattering properties. We consider a condensate in a 
pancake-shaped trap which is flat enough so that the dimension of the cloud in the axial direction is of the order of 
the harmonic oscillator length a z = y/ h/mui z . If the condition a <C a z holds, the system experiences collisions in 
three dimensions and the coupling constant to be used in the 2D NLSE is 

gQ3D =.93d|0o(O)| 2 (4) 

where (j>o(0) = {2-n:a 2 )^ 1 / i is the axial ground-state wave function evaluated at z = 0. This coupling constant has 
been used by a number of authors (see e.g. fl7fl and references therein). When the anisotropy further increases and a 
approaches a z , the collisions start to be influenced by the presence of the trap in the tight z direction and the coupling 
g in such a quasi-2D condensate is to be taken in the form 

gQ2D = 53i)|0o(O)| 2 n - . (5) 

1 + — ^=|ln (2(2ir) 3 / 2 naa z )\ 
a z \J2-K 



This expression was originally derived by Petrov et al. |§,[l8| by studying the scattering function of a system which 
is harmonically confined in the z direction and homogeneous in the {x,y} plane (see also footnote [18] in ||). The 
coupling g in this case depends on density, as is typical of 2D collisions. 

Indeed, when the collisions are in the fully 2D regime (a ^> a z ), the system is described by the coupling 

Anh 2 i 

92D 77-7 — 2TT ' W 

m I In (na ) \ 

as was first derived by Schick for the homogeneous 2D system @. The use of the coupling giD for a system under 
external confinement, involving a dependence on the local density, has been proposed by Shevchenko JlS|] and more 



recently by Kolomeisky et al. |20|. It was rigorously justified by Lieb and coworkers |Ti 

In the following Section we will present results for the condensate wave function ijj(r) as obtained from the NLSE 
(||) using the three couplings gQ3D, and g 2 D- 
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III. EQUILIBRIUM DENSITY PROFILES 



In current experiments on Bose-Einstein condensates of alkali atoms, one of the observables which are most directly 
accessible is the density profile. In particular, it is possible to take in situ images of the cloud, which can be directly 
compared with the results of theoretical models. We give in this Section the predictions for the widths and the shapes 
of the equilibrium density in the {x, y} plane for a condensate in a pancake-shaped trap under different collisional 
regimes as the anisotropy is increased. 

Before proceeding to the numerical solution of the NLSE (Q), it is useful to give a simple analytical estimate of 
the expected width of the cloud. We use for this purpose the Thomas- Fermi approximation, i.e. we neglect the 
kinetic-energy term in Eq. (Q) to obtain the density profile 

\^tf{t)\ 2 = -^(/i - V ext {r))e(ii - V ext (r)). (7) 

Here we have also neglected the spatial dependence of the coupling constant g and used the result of the homogeneous 
system to relate the density to the chemical potential pi] ]. Equation (Q) is valid when the number of atoms in the 
condensate is large and in this regime we have found that its predictions agree well with the full numerical solution 
ofEq. (|). 

The chemical potential in Eq. (Q) is fixed by imposing the normalization condition N = J \ipTF{f)\ 2 d 2 r. By 
inversion of this equation we find an expression for the Thomas- Fermi width of the cloud, Rtf — \/2/i/mw^, in 
terms of the physical parameters of the system. We show in Fig. [I] the width of the cloud as a function of the ratio 
a/a z for a given choice of the number of atoms and of the anisotropy (N — 5 x 10 5 , A = 2 x 10 5 ). We expect the 
quasi-3D model (^) for the coupling to be accurate only for small values for a/a z . It evidently overestimates the size 
of the cloud as a/a z increases. The quasi-2D behavior should be correct for a/a z > 0.1 and is ultimately superseded 
by the purely 2D behavior as a/a z becomes appreciably larger than unity. 

We now turn to the numerical solution of Eq. (||) with the alternative expressions for the coupling. We 

have used the steepest descent method [^2],[l2|, which is known to produce accurate results. A further check of our 
numerical calculations is provided by the virial relation. 

We give illustrations of the ground-state wave function ip(r) predicted by the three models in some relevant cases. 
First of all we have considered the values of particle number, anisotropy parameter, and scattering length as appro- 
priate for 23 Na atoms in the experiment of Gorlitz et al. |Hj (N — 5 x 10 5 , A = 26.33, a — 2.8 nm). In this case 
the system is approaching two-dimensionality for what concerns the confinement (/i ~ 2.08huj z ), but still is 3D for 
collisions (a/a z ~ 3.8 x 1CP 3 ). As is shown in Fig. ||, the quasi-3D and quasi-2D models give almost indistinguishable 
predictions (solid and dotted lines). The fully 2D model would give a quantitatively very different profile (dashed 
line in Fig. @) with a much larger chemical potential (/i = 17.8hu z ), but is evidently inapplicable in this regime of 
parameters. 

For a second case we increase the anisotropy parameter to A = 2 x 10 5 and make the choice a/a z = 0.33 for 
the scattering length. This corresponds to the point of cross-over in the scattering properties from 3D to 2D in a 
condensate where motion in the third direction is completely frozen by the confinement (fi ~ 0.002 huj z ). In this case 
(Fig. ||) the three models predict comparable shapes of the cloud, which is an indication that indeed the condensate 
is entering the regime of 2D collisions. 

Finally, we have further increased the value of the scattering length to a/a z ~ 2.68, keeping the same values for A 
and N. For this choice of parameters the collisions are truly 2D and the fully 2D model should give the most accurate 
prediction for the shape of the cloud, while the other models overstimate its width (Fig. ||) . The calculated chemical 
potential of the gir> model is [i ~ 0.0031 fku z . 

IV. CONCLUSIONS 

In summary, we have considered Bose-Einstein condensates confined inside pancake-shaped traps at zero temper- 
ature within a mean-field description. Depending on the anisotropy of the trap, we have identified three different 
physical regimes for the scattering events: (i) a quasi-3D regime where the axial dimension of the condensate is much 
larger than the scattering length and the collisions are as in a 3D condensate; (ii) a quasi-2D regime where the tight 
harmonic confinement in the z direction begins to influence the scattering events; and (iii) a fully 2D regime where 
collisions are restricted to the {x, y} plane. The atom-atom coupling is different in the three cases and a logarithmic 
dependence on the density arises as the 2D effects start to affect the scattering events. We have adopted a local- 
density approximation to introduce the appropriate atom-atom coupling into the energy functional of the condensate 
and to obtain a non-linear Schrodinger equation for its in-plane wave function. We have given numerical solutions 
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of the NLSE for various choices of parameters corresponding to the three collisional regimes. We have shown that 
the different collisional regimes are reflected in the width of the cloud, which can also be predicted with reasonable 
accuracy by the Thomas-Fermi model if the number of atoms in the condensate is not too low. Our results for the 
fully 2D case are in accord with those reported recently by Lee et al. by a numerical solution of the NLSE with 
the choice <?2Z}(m) f° r the coupling strength. 

In conclusion, the several interesting properties that are expected for flat condensates should be observable when 
not only the condition \l < hoj z is satisfied (the condensate is 2D for what concerns the confinement), but the 
condition a/a z > 1 also holds (the condensate is 2D for what concerns the collisions). The strictly 2D limit can be 
reached in a pancake-shaped Bose- Einstein condensate by making the perpendicular confinement very tight and/or by 
increasing the scattering length, e.g. via Feshbach resonances. In order to understand the significance of interactions 
and high-density effects, high-precision Monte Carlo simulations would be useful as a test of the range of validity 
of the mean-field and local-density approximations. Our analysis should be the starting point for finite-temperature 
calculations, in which to study the role of the interactions with the non-condensed thermal particles and the transitions 
in phase correlations that are expected to occur with increasing temperature. 
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FIG. 1. Transverse width of the condensate (in units of a±_ — sj h/mui±) as a function of a/a z as evaluated in the 
Thomas- Fermi approximation for the three models Q3D (dotted line), Q2D (solid line) and 2D (dashed line). In the inset 
the same curves are plotted over a wide range of the ratio a/a z . The system parameters are N — 5x 10 5 atoms and A = 2 x 10 5 
for the trap anisotropy. 
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FIG. 2. Condensate wave function ip{r) (in units of VN/ax) as a function of r/a±, for a/a z = 3.8 x 10 -3 from the full 
numerical solution of Eq. (Q) for the three models Q3D (dotted line), Q2D (solid line) and 2D (dashed line). The values of the 
particle number and of the anisotropy parameter are indicated in the figure. 
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FIG. 4. The same as in Fig. |[ for a/a z = 2.68 and A = 2 x 10 5 . 



